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Abstract. 

The pointlike impurity scattering is studied for a superconductor with different pairing 
symmetries based on a minimal two band model proposed by S. Raghu et al. [Phys. Rev. B 
77, R220503 (2008)], and the self-consistent Bogoliubov-de Gennes equations. Two intra- 
gap resonance peaks are found at or near the nonmagnetic impurity site for a typical s^^y^- 
symmetry and strong scattering potentials. Such intra-gap resonance peaks are absent for 
the sign unchanged 5- wave and ii^:_,,2-wave order parameter. Thus we propose that a non- 
magnetic impurity with strong scattering potentials can be used to probe the pairing symmetry 
in iron-based superconductors. 
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1. Introduction 

The new family of superconducting materials which contain Fe-As layers has attracted much 
attention since their discovery [IJ. The Fe-As layers are believed to be the conducting 
planes. The band calculations have shown that the Fermi surface includes two hole pockets 
centered at Y points and two electron pockets centered at M points, respectively (21 [31 IHl, 
which is confirmed by angle-resolved photoemission spectroscopy (ARPES) experiments |l5l 
[6l in m 13 [IOl|. A number of experiments have indicated that this family of materials is not 
conventional IfTTTl . Therefore, probing the pairing symmetry is one of the most important 
issues since it can provide us the information of the pairing mechanism. 

So far, the pairing symmetry in iron-based superconductors is still controversial. One 
popular proposal for the pairing symmetry is the 5+-wave gap fill [131 [Ml [IH [161 [13, 
namely, the superconducting gap is extended 5-wave and has opposite sign for the hole and 
electron pockets. Experimentally, the i'-wave gap symmetry is supported by penetration depth 
measurement llTSll . ARPES [[6l [71 [H [91 [ISl and specific heat measurement |[T9l . In particular, 
the gap magnitudes at diff"erent Fermi surface pockets measured by ARPES experiments agree 
well with Sy-iyi = cos cos ky form (here k^ and ky represent the momentum in the unfolded 
Brillouin zone) [[8]|. However, although ARPES experiment is a powerful tool to measure the 
gap magnitude directly, it cannot determine the phase of the superconducting order parameter. 
Detecting the phase of order parameter is an essential and important step to map out the gap 
symmetry. Recently several theoretical works propose different ways to measure the sign 
change of the gap [[201 [ZH [221 [221 [211 • While for multiband materials, probing the phase of 
the order parameter is still a quite challenging task. 

The impurity effect has been an important part in the studies of the superconductiv- 
ity [[25l . The effect of the impurity scattering is sensitive on the symmetry of the order pa- 
rameter, thus it is a useful tool to probe the pairing symmetry. A well-known result of i^-wave 
pairing in the cuprates is the zero bias peak in the local density of states (LDOS) near the non- 
magnetic impurity site. For the case of iron-based materials, it was proposed that the impurity 
scattering should be sensitive to the sign change of the 5+-symmetry order parameter [FTIl. 
Furthermore, the averaged effect of an ensemble of impurities on the thermodynamic proper- 
ties has been studied in iron-based materials [[26l[27l[28l[29l[30l . The impurity-induced intra- 
gap state forms due to the inter-band scattering [[29ll for the sign-reversing fully gapped order 
parameter symmetry. We note that the calculation in Ref [ 291 is based on a phenomenological 
approach, namely, they treat the intra-band scattering and inter-band scattering intensity as 
input parameters and consider an impurity concentration. However, such intra-gap state is not 
observed in scanning tunneling microscopy (STM) experiments [[3Tl[32l[33l[34l . In fact, the 
impurity-induced intra-gap states are expected to exist only at or near the impurity site. To 
study the possibility and property of such states one effective method is to put one or several 
pointlike impurities in real space and calculate the LDOS at or near the impurity site. Basi- 
cally two methods can be used, namely, the analytic T-matrix approach and the self-consistent 
Bogoliubov-de Gennes (BdG) equations. Based on the first method and taking into account 
the s+ symmetry phenomenologically, very recently Ref. |[35l calculated the LDOS near the 
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impurity site and do find the sharp resonance peaks inside the gap. The peaks will disappear 
if the phase of the order parameter is taken to be the same. The intra-gap bound states are 
actually located at the positive and negative energies and the LDOS tends to be zero at low 
energies. The impurity effect is expected to be sensitive to the detailed band structure [|25l . 
The sensitivity maybe more remarkable for the multi-band materials. Since the calculation in 
Ref. BSil is based on a band structure which artificially fits the ARPES experiments [|36l[37ll . 
in the present paper we wish to examine the LDOS at and near the impurity site based on 
a different two-band model as proposed by Raghu et al. If38]| . and the self-consistent BdG 
equations. The advantage of the BdG technique is that it calculates the order parameter in 
a fully self-consistent way. The gap magnitude is expected to be suppressed at and near the 
impurity site. It is proposed that such suppression effect needs to be taken into account to 
study the local impurity states ll39l . i.e., the suppression effect would enhance the intensity of 
the bound states and the peak width and position also depend on the amplitude of the suppres- 
sion. For the case of the s^iyi pairing symmetry, our numerical results show that the intra-gap 
states exist as long as the scattering potential is strong enough, but the detailed spectra of 
the local impurity states are quite different and sensitive to the parameters for obtaining the 
band structure. For a nonmagnetic impurity with positive scattering potential, we obtain two 
intra-gap resonance peaks. The resonance peaks are located near the coherent peaks for weak 
potentials and move to lower energy as the potential increases. The weight of the resonance 
peak with positive energy becomes very pronounced as the scattering potential reaching the 
unitary limit. For the cases of negative scattering potentials, the intra-gap resonance peaks are 
moving away from the coherent peaks and are closer to the Fermi energy. The position of the 
resonance peaks are robust and depend weakly on the scattering potentials. All other features 
remain practically the same. We also studied the LDOS spectra for the cases of the sign un- 
changed i'-wave and nodal d^2_y2-wa.\e and no intra-gap peaks are observed. Experimentally 
the LDOS can be measured by STM experiments, and thus is a useful tool to probe the pairing 
symmetry of this family of compounds. 

The paper is organized as follows. In Sec. II, we introduce the model and work out the 
formalism. In Sec. Ill, we persent the numerical calculations and discuss the obtained results. 
Finally, we give a brief summary in Sec. IV. 

2. Model and formalism 

As we mentioned above, iron-based superconducting materials have a layered structure. In the 
FeAs layer there are two Fe ions per unit cell by taking into account As ions located above or 
below the Fe-Fe plane. Currently no consensus has been reached on what is the suitable model 
to capture the essential physics [IM |40l SB |42l |43l HH SSI SH. Some groups argue all five 
3d orbitals have to be considered to construct a minimal model [ l4ni42| . On the other hand, 
the band calculations have also shown that the main features of the bands that determine the 
Fermi surface are the two orbitals and dy^ [|47ll48l . Thus the two band model is proposed 
by several groups [l38l l45l l46l . Because we are concerned with the low energy physics in 
the present work, we expect that the qualitative result does not depend on the detailed band 
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Figure 1. The Fermi surface in the extended Brillouin zone and folded Brillouin zone 
respectively. The black and green dashed lines are nodal lines for .s^2,,2-wave, and dx^-y-^- 
wave pairing symmetry, respectively. The red and blue colors represent the hole and electron 
pockets, respectively. 



Structure. We here use the two band model considering the two degenerate orbitals dxz and dyZ 
per site Il38ll . On the two-dimensional square lattice, the model Hamiltonian reads, 



H - Ht + H/^ + Hjfnp, 
where is the hopping term. 



(1) 



' i/io- • (2) 

i/i }V(T lf.l<T 

where /, j are the site indices and yu, v = 1, 2 are the orbital indices, to is the chemical potential. 
//a is the pairing term. 



Ha= J]{Ai,jycl^cl^ + h.c.). 



(3) 



Himp is the impurity term, which is expressed by. 



imp 



(4) 



where i^ represents for the impurity site and V^^y is the scattering potential. 

The Hamiltonian can be diagonalized by solving the BdG equations self-consistently. 
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where the Hamiltonian //,-u,v is expressed by. 



(5) 



(6) 



The superconducting order parameter and the local electron density n,^ satisfy the 
following self-consistent conditions. 



2KbT 



(7) 
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n n 

Here Vt^jv is the pairing strength and f{x) is the Fermi distribution function. 
The LDOS is expressed by, 

P,M = Y}\ulJ^6{En -oj) + \vl^?6{En + oj)], (9) 

where the delta function 6(x) is taken as Y/7t{x^ + Y^), with Y = 0.004. The supercell technical 
is used to calculated the LDOS. 

In the following calculation, we use the hopping constant suggested by Ref. If38]| . where 
?i_3 are hoping constants between the same orbital, expressed by. 



tnj±n - ti2,i±y2 - ti - —l.Q, (10) 

ti2,i±x2 = ^il,(±.vl = h = 1-3, (11) 

ti;,J±,±y, = t3 = -0.S5 (ju=h2). (12) 
?4 is the inter-orbital hopping constant, expressed by, 

ti,i,i±(x+y)v = t4= 0.85 (yU ^ v), (13) 

tif,j±(x-y)v = -U = -0.85 (m + V). (14) 



?o is determined by the doping density. Throughout the work, the energy is measured in units 
oi\t\\ = 1. The Fermi surface with filling electron density 2.1 per site is plotted in Fig.[IJa). 
Alternatively, if we consider one unit cell containing two Fe ions, the Fermi surface in the 
reduced Brillouin zone is plotted in Fig. [2b). 

The pairing symmetry is determined by the pairing potential Vi^jv We have carried 
out extensive calculations to search for favorable pairing symmetry. To summarize, the 
on-site potential, the nearest-neighbor (NN) potential, and the next-nearest-neighbor (NNN) 
potential, will produce isotropic j'-wave (A = Aq), J_^2_y2-wave [A = 2Ao(cos - cos ky)], and 
s^2y2-vja\e (A = 4Ao cos k^ cos ky) pairing symmetry, respectively. In the present work, we will 
focus on the characteristic of s+ symmetry. We consider the pairing between the same orbital 
of the NNN site. Without the impurity, self-consistent calculation verified that this kind of 
pairing produces 5-wave gap symmetry, i.e., A,^ = Aq. Transforming the Hamiltonian 
to the momentum space and diagonalizing the 4x4 Hamiltonian, it is easy to verify that the 
pairing symmetry is exactly the same with the 5 v2y2 -symmetry proposed by Ref. (49], namely, 
the gap function in the extended Brillouin zone has cos k^ cos ky form, and has the form of 
cos k'^ + cos ky in the reduced Brillouin zone. The nodal line of the gap function is shown in 
Fig. 1. 

In the following presented results, we choose the filling electron density n = 2.1 per site 
(electron doped samples with doping 6 = 0.1) and pairing potential V = 1.2. We consider 
the case for the intra-orbital scattering, namely, V^^y = KcJ^jv The parameters are chosen 
just for illustration. We have checked numerically that our main results are not sensitive to 
the parameters. The numerical calculation is performed on 20 x 20 lattice with the periodic 
boundary conditions. The impurity is put at the site (10, 10). A 80 x 80 supercell is taken to 
calculate the LDOS. 
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Figure 2. Amplitudes of the order parameter in presence of a single nonmagnetic impurity 
with y, - +5. 

3. Results and discussion 

We first illustrate the feature of single nonmagnetic impurity scattering with Sx2y2 pairing 
symmetry. The scattering potential Vs can be estimated roughly by comparing the energy 
of the impurity atoms and that of Fe atoms. In real materials, both negative and positive 
scattering potentials are possible. Thus here we discuss both cases. As it is unclear what an 
accurate value of the potential strength is, we consider a range of potential values. 

We plot the self-consistently determined order parameter amplitudes on the lattice 
(A,- = |liy^A;^j^) in Fig. |2] for = ±5. As seen, the order parameter is suppressed near 
the impurity and quite small at the impurity site. It will recover to the uniform value at about 
2 ~ 3 lattice spacings from the impurity site. Note that the suppression effect is more notable 
for the negative potential. 

The LDOS spectra for different scattering potentials at and near the impurity sites are 
plotted in Fig. [31 The left panels are the negative potential cases. Obvious impurity states 
(denoted by the red arrows) exist at the energy about ±0.57A (here A is the energy of the 
superconducting coherent peak, denoted by the green arrows). The energy of the impurity 
resonance peaks depends weakly on the scattering potential V^. At the impurity site, for 
small scattering potential V^, two sharp resonances peaks inside the gap are shown, and the 
resonance peaks are suppressed as Vs increases and the LDOS tends to zero at the unitary 
limit. On the other hand, near the impurity site, the intra-gap resonance peaks show up and 
the intensity increases as the potential increases. The resonance peaks are remarkable at the 
NNN site. In the unitary limit, a very sharp and strong resonance peak appears at the positive 
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Figure 3. The LDOS spectra in presence of a single nonmagnetic impurity for different 
scattering potentials. The black lines are the LDOS at the impurity site. The red and blue lines 
are the spectra at NN and NNN sites to the impurity. The left and right panels are for negative 
and positive scattering potential, respectively. The green line in panel (a) is the density of 
states spectrum without the impurity. The green and red arrows indicate the superconducting 
coherent peak and intra-gap resonance peaks, respectively. 

energy for the LDOS at the NNN site of the impurity site. This situation may qualitatively 
correspond to the case for a charged Ba^^ impurity with a negative or attractive potential seen 
by electrons on the FeAs layer. This type of impurity may be present in the FeAs layer when 
BaFe2-vCOvAs2 is used for STM experiments [|50l . 

The right panels show the effect of a positive scattering potential impurity. We do not see 
clear resonance peaks at the impurity site for the weak potential, which is different from that 
of the negative potential. Near the impurity site, similar to the case of the negative potential. 
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Figure 4. The LDOS spectra in presence of a single nonmagnetic impurity for on-site i-wave 
symmetry and d^-^-x^- symmetry with Vj - +5, respectively. The black lines are the LDOS at 
the impurity site. The red and blue lines are the spectra at NN and NNN sites to the impurity. 
The green lines are density of states spectra without the impurity. 



the intra-gap resonance peaks show up as the potential increases, while the intensities seem 
to be smaller than those of the negative potential impurities. The resonance peaks are more 
remarkable at the NNN site to the impurity, which is similar to that of the negative potential. 
We also point out that the positions of the two resonance peaks are much closer to the coherent 
peaks for weak scattering potentials. The intensity of the resonance peaks with positive energy 
increases and the positions move away from the superconducting coherent peaks to lower 
energies as the scattering potential increases. For a quite strong impurity scattering potential 
{Vs = 10), the intra-gap peaks are located at the energy about +0.7A. The spectrum is quite 
similar to that of the negative potential (Vg = -10). 

On the whole, we present the LDOS spectra taking into account the impurity scattering 
for 5;f2j,2 -pairing symmetry. The intra-gap states for 5;(.2j.2 -pairing symmetry obtained here are 
quite robust for a strong impurity. Experimentally, the intra-gap bound states may be detected 
by comparing the LDOS spectrum far away the impurity site with that near the impurity site. 
The intra-gap peaks will show up near the impurity site. We also verified numerically that the 
main results will not change for many impurities or the impurity scattering potential is not 
only purely on-site but also extends with a few lattice sites. While for quite weak scattering 
potentials, the intensity of the intra-gap peaks would decrease and disappear as the potentials 
are small enough. We have checked that no obvious intra-gap features could be observed as 

We now discuss the impurity scattering effect for different pairing symmetries. The 
LDOS spectra at and near the impurity site for the sign unchanged 5- wave (on-site pairing) 
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and (i_^2_y2-wave (NN pairing) with = ±5 are shown in Fig. HI For the case of the sign 
unchanged i'-wave pairing symmetry, we can see clearly two superconducting coherent peaks 
with the energies equal to those of the density of states without the impurity. The intensity 
of the coherent peaks are suppressed at the impurity site and its NN site. No intra-gap peaks 
are observed. For the case of the J-wave pairing symmetry, two intra-gap plateaus can be 
observed inside the gap. The spectra are significantly different with those of the cuprates. 
In fact, here the calculations are based on the two-band model. For the J^2_j,2-wave order 
parameter, the hole pockets around the F point are nodal and the electron pockets around 
the M point are nodeless. As a result, the low energy spectra are mainly contributed by the 
hole pockets, which lead to the intra-gap plateaus. The interband scattering may enhance the 
intensity of the plateaus while no intra-gap bound states are observed for d^2_y2 symmetry. 

On the whole, we present the LDOS spectra taking into account the impurity scattering 
for different pairing symmetries. The intra-gap states for 5^-2^,2 -pairing symmetry obtained 
here are quite robust for a strong impurity. The electron operator for each orbital is a 
linear combination of quasiparticle operators from all bands, therefore the four Fermi surface 
pockets in Fig. 1(b) are contributed by the hybridigation of the two orbitals. With this 
observation, one can find out that even the intra-orbital impurity scattering will give rise to the 
scattering between different bands with the scattering strength weighted by the Bogoliubov 
amplitude IfSTI . Comparing the LDOS spectra for the s^iyi pairing symmetry shown in Fig.3 
and those for the sign-unchanged 5-wave symmetry shown in Figs. 4(a) and 4(b), we expect 
that the intra-gap bound state for the s^iyi pairing symmetry is due to sign change of the order 
parameter between the different Fermi surface pockets. As a result, we propose that this can 
be used to detect the sign reversal of the order parameter. However, we stress that the sign 
change of the order parameter is just one of the essential conditions for the bound states but 
not sufficient. The scattering is weighted by the Bogoliubov amplitude so the band structure 
and the detailed scattering potential form is also important. For the current model, the d- 
wave pairing symmetry will not produce the bound states, as seen in Figs. 4(c) and 4(d). We 
here focus our discussion on the properties of the s+ pairing symmetry. The reason for the 
absence of the bound states for J- wave pairing symmetry can be investigated further through 
the momentum space scattering potential based on the T-matrix method while this issue is not 
concerned here. 

Comparing all the relevant results with different pairing symmetries, we propose that one 
strong nonmagnetic impurity can be used to probe the pairing symmetry and detect the sign 
reversal of the order parameter. 

4. Summary 

In summary, we have studied the point-like single non-magnetic impurity scattering effect in 
iron-based superconductors based on self-consistent BdG equations. The LDOS at and near 
the impurity site are calculated and the spectra have their unique features for different pairing 
symmetry. Especially, for the case of the ^•^-2,,2-wave symmetry, two intra-gap peaks show up 
as the impurity scattering potential is strong enough. This may also be easily detected by the 
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STM experiments through substituting the Fe atom with a nonmagnetic atoms (e.g., Zn or Ba 
atoms etc.). As a result, we propose that the non-magnetic impurity can be used to probe the 
pairing symmetry of iron based superconductors. 

Note added- After we posted our work in Arxive, we note that another group also did 
similar calculations [53] based on BdG equations and T-matrix method using the same two 
band model proposed by S. Raghu et al. Their results are similar to ours. 
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